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Abstract. We study the nonequilibrium Casimir-Polder force on an atom prepared 
in an incoherent superposition of internal energy-eigenstates, which is placed in a 
magnetoelectric environment of nonuniform temperature. After solving the coupled 
atom-field dynamics within the framework of macroscopic quantum electrodynamics, 
we derive a general expression for the thermal Casimir-Polder force. 



PACS numbers: 12.20.-m, 42.50.Ct, 42.50.Nn, 
1. Introduction 

The thermal fluctuations of the electromagnetic fleld present at flnite temperature may 
interact with a single atom or molecule, resulting in the thermal Casimir-Polder (CP) 
force. It has been studied theoretically via linear response theory [IJ, Lifshitz theory 
[2], quantum electrodynamics (QED) [3] as well as heuristic generalisations of the zero- 
temperature result |1] and experimentally by spectroscopic means [5]. 

Thermal CP forces in nonequilibrium scenarios have recently received an increased 
attention, where two cases can be distinguished. A nonequilibrium environment can 
be realised when an atom interacts with a body whose temperature differs from the 
environment temperature |6]. Alternatively, situations have been considered where the 
environment is at global thermal equilibrium with a uniform temperature, but one or two 
present atoms are not fully thermalised with their local environment [3, IE]- In the present 
paper, we generalise these two special cases to the full nonequilibrium situation of an 
atom in a nonequilibrium state placed within an arbitrary magnetoelectric environment 
of nonuniform temperature. We use macroscopic QED (recalled in Sec. [2]) to calculate 
the CP force from the thermal average of the quantum Lorentz force (Sec. Ej) followed 
by a short summary (Sec. H]). 



2. Macroscopic quantum electrodynamics at finite temperature 

Consider an environment of dispersive and absorbing magnetoelectrics of (relative) elec- 
tric permittivity e{r,u}) and a magnetic permeability fi{r,u), which are both satisfying 
the Kramers-Kronig relations. This environment together with its electromagnetic fleld 
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can be characterised by a Hamiltonian Hp = J2x=e m I Io° ^ /a(^' ^) ' /a(^) ^) 
P [in], where the bosonic variables fx and /| are associated with the noise polarisation 
(A = e) and magnetisation (A = m) of the media. When a single atom or molecule 
with Hamiltonian = ii^n I n)(n I {En- eigenenergies, |n): molecular eigenstates) is 
placed at position in this environment, it interacts with the electromagnetic field via 
an electric-dipole Hamiltonian Haf = - Y.m,n dmn- E{rA)Aran {dmn = {m\d\n), Amn = 
\m){n\), so that the total Hamiltonian of the system reads H = Ha + Hp + Haf- 

The electric field can be expressed in terms of the dynamical variables according to 

Eir)= / duE{r,u) + 'R.c.= du d^r Gxir,r' ,u)- fxir' ,u) + R.c, (1) 

with the coefficients Ga being related to the classical Green tensor, G, by 



I h 

Ge{r,r',uj) =i^\ lrae{r' ,uj) G{r,r' ,uj), (2) 



V TT^o 




G„.(...'..) ^J^i^[V'xG(/,.,.)] . (3) 



The Green tensor is the unique solution to the Helmholtz equation 



2 



V X — — -V X e(r,uj) G(r,r',uj) = S(r-r') (4) 

with G(r, r', w) — )■ for |r — r'| — > oo, where the above definitions imply 

J2 I dhG,{r,s,u)-Gl\r',s,co) = ^ uHmGir,r' ,co). (5) 

A=e,m 

For a nonuniform temperature T = T(r) the thermal state of the environment 
is described by pr = exp{-/ d^r ^FW/[fcBT(r')]}/ Tr(exp{- / d^r HF{r)/[kBT{r)]}) 
[Hp = J d^rHp{r), k^'- Boltzmann constant]. The relevant nonvanishing thermal 
averages of the dynamical variables are thus given by 

{flir, u)Mr', co')) = nrir, u)5xx'Sir - r')5(a; - to'), (6) 
{fx{r,io)fi,{r',uj')) = [nT{r,io) + l]6xx'S{r-r')6{io-uj') (7) 

where nT{r,u}) = l/{exp{f}Lo/[kBT{r)]} — 1) is the average thermal photon number. 



3. The thermal Casimir— Polder force 



The thermal CP force on an atom prepared in an incoherent superposition of internal 
energy-eigenstates can in electric dipole approximation be found from the average 
Lorentz force [101 E] -^(^a, i) = { [^d-E{r)~\^^^^y In order to evaluate this expression, 
one needs to solve the coupled atom-field dynamics. Using the Hamiltonian given in 
Sec. [21 one finds 

-Amn = ^^mnAmn + ^ ^ {dnkA-mk ~ '^fcm^fcn) 'EiVA)-, (8) 

k 
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fx{r,u}) = -iufxir,u) + ^^dmn-Gl{rA,r,u)Amn- (9) 

m,n 

We eliminate the field by formally integrating the second of these equations and 
substituting the result back into the first one, which we arrange in normal ordering. 
After invoking the integral relation (|5]), one obtains 

k 

+e-''^*[i™fc(t)d„fc - Akn{t)dkm] -MrA.uj)] + Z^n{t), (10) 

with 

Z^nit) = yY. r diouj'dM-lmG{rA,rA,io)- T dr{e-^-(*-^) 

^ k,l,j ° ° 

-d^jA^j{t)]AM{T) + e^"(*-")ifez(r) [d^jA^jit) - dj^Anjit)] } (11) 

denoting the zero-point contribution to the internal atomic dynamics. The thermal 
contribution can be determined iteratively by substituting the self-consistent solution 



I. ^ 



X [Amk{T)dnk - Akn{T)dkrr,] ■ [E{rA, w)e-^^" + H.C.] (12) 

to the truncated Eq. ffTOj) [without Zmn{t)] back into Eq. flTOl) and taking thermal 
expectation values with the aid of Eqs. (|H]) and ([7]). This leads to a closed system 
of equations {Amn)= i^jJmn{Amn) + {Zmn)+{Tmn) with a thermal contribution 

(f„„(t)>=-^V Hdu V [ d'rnT{r,co) T dr [e-'-(*--) + e^-(*--)] (13) 

^ k,l X=e,m'^ 

x{e''^-'^(*--)d„,-GA(r^,r,cu)-Gf(r^,r,a;)-[dH(i™K^)>-^i^rn(Afe(r)>] 
_e-.n(*-)d,^.G;,(rA,r,c^)-Gf(rA,r,c^)-[d^(^(^)>-^i^fc(^^n(r)>]}. 

Assuming the atom-field coupling to be sufficiently weak, we can apply the Markov 
approximation by writing <^Amn(T)) — e^"^'""*^*^'^^(y4mn(^)) and letting the lower limit of 
the time integrals tend to minus infinity, so that dr e''^'^*^^-' ~ 'k5{x) + iP / x. For a 
nondegenerate system, the off-diagonal elements of the (internal) atomic density matrix 
(T decouple from each other as well as from the diagonal ones, and one finds that internal 
atomic dynamics follows the rate equations 

(^nn{t) = - r„a„„(t) + ^ TknCTkkit), (14) 

k 

arnn{t)=[-iUrnn - ^(Xrn + 'rn)/2]arnn(t) loi 111 ^ U (15) 

(c"mn = {m\a\n) = (^Anm))- As follows from Eqs. (ITTj) and (fT3!) . the total loss rates read 
Tn = ^r„fe = ^ ^a;^^^e(a;„fc)d„fc-ImG(rA,n4, \unk\)-dkn 

k k 
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k A=e,m 

xdnk-GxirA,r, \u)nk\)-G*^{rA,r, \oJnk\)-dkn (16) 
and the shifts of atomic transition frequencies Umn = ^mn + S^m — S^n are given by 



du CO 



X=e,m 



^nk — ^ 

nT{r,uj) ^ nT{r,uj) 

^nk — ^ ^nk + 



xdnk-Gx{rA,r,uj)-G*^{rA,r,uj)-dkn. (17) 

Note that we have replaced the Green tensor with its scattering part G*-^-* for the zero- 
point frequency shift since the free-space zero-point Lamb shift is thought to be aheady 
included in the bare transition frequencies Umn- 

Having solved the atom-field dynamics, we can now evaluate the average Lorentz 
force by using Eqs. ([1]), ([6]), ([7]), f|T2|) and the solution to Eq. ([9]). After invoking the 
correlation function (v4„„(t)v4fc/(r)) = 5„fc(^mi(^))e'^'""^*~^-' [^^mn = '^m„ + i(r„-Fr„)/2] 
which follows from Eq. f|T5|) via the quantum regression theorem [12], one obtains 



^(^A,t) = -^^/ I d^sVdnk-Gx{r,s,uj)-G*^{rA,s,uj)-dkn 



n,k 



X=e,m 



X 



[ dr(i„„(r)){nT( 

^0 



s, a; e 



i(w+r2„fc)(t-r) 



+ [nT{s,u)) + l]e 



-i(tj-n„fe)(t-T) 



r=rA 



} +C.C. 



Using the Markov approximation in the form dr (A„„(r)). . . ~ (A„„()f:)) dr . . ., 
the thermal CP force for an atom prepared in an incoherent superposition of internal 
energy-eigenstates is given by i^(r4, t) = (T„„(t)i^„(rA) with force components 



(19) 



^ — ^nk ^ + ^nfc 

X Vd„fc-GA(r, s,uj)-G*^{rA, s,u)-dkn\^^^^ + c.c. 

An important special case is that of an equilibrium environment at temperature T. 
Applying Eq. ^ to the s-integral in Eq. f[T^ and transforming the tu-integral by means 
of contour- integral techniques, one obtains in the perturbative limit (VLnk — ^nk) 



FnivA) = -fioknTj^i^ - -,5m)^l'^ATr[arriiO-G^'\rA,rA,i^N)] (20) 

N=0 

+/^oX^{0(Wnfc)h(w„fc) + 1] - Q{uJkn)n{uJkn)}ujlk^Adnk-'R'eG^^\rA,rA,^nk)-dkn 



{C,N = ^Trk^TN/h, Matsubara frequencies), where the atomic polarisability reads 

d-nkdkn , dknd^k 



OLniUJ 



lim — 



-Unk -OJ -le 



+ 



-ujnk + w + ie 



(21) 
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This is in agreement with previous results |7]. Equations (IT^-( |T6ll show that in the 
long-time hmit one haso-(t — )■ oo) = (3-r = exp[— _f/'yi/(fcB^)]/Trexp[— ^'^/(/cbT)], and the 
Lifshitz force F(r^,t ^ oo) = -^^okBTY.N=o{^ - ^M^^^OOVa Tr G(^)(r^, r4, i^iv) 
is recovered with a thermal polarisability aT(i^) = crr,nn«n(i^)- 

4. Summary 

We have used macroscopic QED to obtain a general expression for the CP force on 
an atom which is placed in an arbitrary magnetoelectric environment of nonuniform 
temperature and whose initial internal state may be an arbitrary superposition of 
internal energy-eigenstates. Our result reduces to previous ones in the special case 
of uniform temperature. 
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